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A code for predicting supersonic jet broadband shock-associated noise was assessed us- 
ing a database containing noise measurements of a jet issuing from a convergent nozzle. 

The jet was operated at 24 conditions covering six fully expanded Mach numbers with 
four total temperature ratios. To enable comparisons of the predicted shock-associated 
noise component spectra with data, the measured total jet noise spectra were separated 
into mixing noise and shock-associated noise component spectra. Comparisons between 
predicted and measured shock-associated noise component spectra were used to identify 
deficiencies in the prediction model. Proposed revisions to the model, based on a study 
of the overall sound pressure levels for the shock-associated noise component of the mea- 
sured data, a sensitivity analysis of the model parameters with emphasis on the definition 
of the convection velocity parameter, and a least-squares fit of the predicted to the mea- 
sured shock-associated noise component spectra, resulted in a new definition for the source 
strength spectrum in the model. An error analysis showed that the average error in the 
predicted spectra was reduced by as much as 3.5 dB for the revised model relative to the 
average error for the original model. 

I. Introduction 

A jet noise prediction model must compute the radiated noise from all the sources that exist throughout 
the range of useful jet operating conditions. At lower speeds, the noise may be characterized as coming 
from the fine-scale turbulence in the flow. As the speed increases, the noise due to coherent, large-scale 
turbulence becomes significant starting in the directions off the downstream jet axis. Near sonic conditions 
and above, pressure imbalances in the flow create shocks. The interaction between turbulence and shocks 
creates an additional source of noise that radiates over a broad range of higher frequencies and predominantly 
in the direction upstream of the nozzle exit. If coherent feedback occurs, then intense screeching can also be 
present. These are all potential sources of radiated jet noise. An approach to making jet noise predictions 
is to model the sources separately. Then assuming incoherent sources, the radiated noise from each source 
is combined in the far field to obtain the total radiated noise. Once a modeling tool has been developed, it 
must be validated by comparing predictions to experimentally measured data. 

As part of a larger study to assess NASA’s capability to predict aircraft noise, empirical and acoustic 
analogy based tools were used to predict jet noise and the results were compared to measured data over a wide 
range of operating conditions. 1 The jet experiments used nozzle geometries including convergent, convergent- 
divergent, internal mixer, separate flow, and chevron configurations operating with Mach numbers from 0.3 
to 2.0 and temperatures from 0.8 to 3.0 times the ambient temperatures. The predicted jet noise included 
both mixing noise, for all jet operating conditions, and broadband shock-associated noise, for supersonic jet 
operating conditions. The general conclusions were that the empirical tools could make successful overall 
sound pressure level predictions for a broad range of nozzle geometries and operating conditions while the 
acoustic analogy tools were successful for a very limited range of operating conditions using a single stream 
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nozzle. No attempt was made in the above study to separately assess the ability to predict the broadband 
shock-associated noise generated during supersonic jet flow conditions. 

Dahl 2 recently reported separate assessments of jet mixing noise and broadband shock-associated noise 
prediction codes. Included in that report was an assessment of the SAE broadband shock-associated noise 
prediction model that is documented in SAE ARP 876E 3 and is also included within NASA’s Aircraft 
Noise Prediction Program (ANOPP). 4 The measured jet noise data used in the assessment were total far- 
held noise spectra decomposed into component spectra for mixing noise and shock-associated noise. The 
predicted shock-associated noise spectra were compared to the shock-associated noise component of the 
measured spectra. Favorable agreements were noted for the prediction codes with the exception of the 
SAE model code. This code gave predicted spectra that compared favorably to the total noise spectra 
in the upstream direction and became increasingly inaccurate at higher jet temperatures and downstream 
directions. According to the SAE document, the total predicted jet noise is the sum of the predictions from 
the mixing noise model and the shock-associated noise model. In all cases shown in Dahl, 2 the addition 
of mixing noise to the SAE shock-associated noise code output would lead to inaccurate results especially 
at lower frequencies where mixing noise dominates the spectrum. This paper examines the SAE shock- 
associated noise prediction model using the shock-associated noise component of the measured data and 
proposes changes to the model. The model would then meet the intent stated in the SAE document of only 
predicting the shock-associated noise component of supersonic jet noise and that the total predicted radiated 
noise would be obtained by adding the predictions from a jet mixing noise model. 


II. Experimental Data 


The experimental jet noise spectral data used in this study were acquired during an extensive test program 
conducted using the Small Hot Jet Acoustic Rig (SHJAR) at the NASA Glenn Research Center. The SHJAR 
is a single-stream, hot jet rig installed within the Aero- Acoustic Propulsion Laboratory, a 65-foot geodesic 
dome lined with acoustic absorbing material to provide anechoic conditions for noise testing at all frequencies 
above 200 Hz. More information on this facility may be found in Bridges & Brown 5 and other reports describe 
the types of data and the quality of the data collected during testing in this facility. 6-8 

Khavaran & Bridges 9 discuss details of the test results from the program to measure the noise from jets 
at a variety of nozzle pressure ratios and total temperatures. Tables are given that describe the different 
nozzles used and the variety of operating conditions tested. A subset of that data is used here consisting of 
noise data from supersonic jets issuing from a 2.0-inch exit diameter, circular, convergent nozzle. The nozzle, 
designated as smc021, is a modified version of the smcOOO baseline convergent nozzle with a 5-degree conic 
contraction and a 0.04-inch thick lip, see Figure 1. The smc021 nozzle had a sharp edge lip with notches cut 
into the edge. These nozzles are shown in Figure 2. The purpose of the modifications was to minimize the 
screech noise observed in the smcOOO nozzle spectral data. 10 




Figure 2. Nozzles: smcOOO on the left, smc021 on the 
right. 


To assess shock-containing jet noise prediction models, the data were collected using the smc021 nozzle 
operated at the 24 conditions listed in Table 1. The operating conditions are grouped according to constant 
total temperatures. However, the operating condition parameters were recorded at the time the acoustic 
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Table 1. Test conditions for a convergent nozzle designated smc021 grouped according to total temperature. 


Rdg 

Ti/Too 

NPR 

Mj 

V j / C 0 o 

P 

7 3 

Too(K) 

Poo (Pa) 

Coo (nr/s) 

2079 

1.01 

2.005 

1.049 

0.954 

0.315 

1.4000 

279.4 

97789 

335.04 

2080 

1.01 

2.358 

1.178 

1.047 

0.623 

1.4000 

279.2 

97783 

334.92 

2081 

1.01 

3.147 

1.392 

1.186 

0.969 

1.4000 

279.2 

97788 

334.93 

2082 

1.01 

3.636 

1.493 

1.245 

1.109 

1.4000 

279.2 

97790 

334.93 

2083 

1.01 

4.647 

1.660 

1.338 

1.325 

1.4000 

279.2 

97796 

334.89 

2085 

1.01 

5.708 

1.796 

1.406 

1.491 

1.4000 

279.1 

97794 

334.87 

2088 

1.77 

2.004 

1.050 

1.264 

0.319 

1.3951 

280.1 

97814 

335.44 

2089 

1.77 

2.365 

1.182 

1.388 

0.630 

1.3963 

280.2 

97817 

335.52 

2090 

1.76 

3.143 

1.392 

1.568 

0.968 

1.3978 

280.2 

97826 

335.49 

2091 

1.76 

3.627 

1.492 

1.647 

1.108 

1.3983 

280.2 

97831 

335.50 

2092 

1.77 

4.644 

1.660 

1.775 

1.325 

1.3990 

280.0 

97826 

335.37 

2093 

1.76 

5.674 

1.792 

1.858 

1.487 

1.3994 

279.8 

97814 

335.30 

2100 

2.16 

1.989 

1.045 

1.390 

0.305 

1.3872 

279.9 

97851 

335.34 

2098 

2.16 

2.356 

1.181 

1.533 

0.628 

1.3896 

279.9 

97842 

335.34 

2097 

2.16 

3.130 

1.390 

1.734 

0.966 

1.3929 

279.8 

97837 

335.25 

2096 

2.16 

3.629 

1.494 

1.828 

1.109 

1.3942 

279.6 

97839 

335.18 

2095 

2.16 

4.638 

1.659 

1.959 

1.324 

1.3961 

279.6 

97825 

335.16 

2094 

2.16 

5.685 

1.794 

2.056 

1.489 

1.3973 

279.7 

97830 

335.23 

2102 

2.66 

1.987 

1.048 

1.542 

0.314 

1.3736 

280.2 

97865 

335.54 

2103 

2.65 

2.327 

1.175 

1.691 

0.617 

1.3769 

280.3 

97857 

335.58 

2104 

2.66 

3.120 

1.391 

1.926 

0.967 

1.3823 

280.3 

97859 

335.59 

2105 

2.65 

3.597 

1.490 

2.020 

1.105 

1.3849 

280.4 

97864 

335.61 

2106 

2.64 

4.630 

1.660 

2.168 

1.325 

1.3889 

280.4 

97857 

335.61 

2107 

2.65 

5.638 

1.789 

2.276 

1.484 

1.3913 

280.3 

97854 

335.58 


data were measured, hence, the numbers in the table reflect the slight variations in conditions that occurred 
during the testing. The table column headings are given as follows: 


Rdg - designates the experimental reading number used to label the various operating conditions. 
Tt/Too - the nozzle flow total temperature to ambient temperature ratio. 

NPR - the nozzle pressure ratio defined as the nozzle flow total pressure divided by the ambient 
pressure. 

Mj - the fully expanded jet Mach number computed using the equation 


Mj = 



( 1 ) 


Vj/c oo - the acoustic Mach number is the ratio of the fully expanded jet velocity to the ambient 
speed of sound. The fully expanded jet velocity is computed from Mj using the fully expanded 
speed of sound 

c i = \/ R Tj ( 2 ) 

and the fully expanded temperature 


Tj = 


1 + 


T t 


( 3 ) 
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7 j - the ratio of specific heats at the fully expanded jet conditions. It is computed following the 
procedure given in the Appendix A. 

/3 - the shock noise parameter 

P = v/M? ^ 1 (4) 

Tqo - the ambient temperature in degrees Kelvin. 

Poo - the ambient pressure in Pascals. 

Coo - the ambient speed of sound in meters per second computed from the equation 

Coo = ^JlRToo (5) 

where 7 = 1.4 and R = 286.959 m 2 /(s 2 K). 

III. Description of Measured Data and Analysis 

The measurements of the noise radiated from the jet were made by 24 microphones. These microphones 
were located on an arc 50 nozzle exit diameters from the center of the nozzle exit area and placed at 24 
equally-spaced angles. The angles were from 50 to 165 degrees as measured from the upstream inlet centerline 
to the nozzle. 

As described in Bridges et al . 11 and in Khavaran & Bridges , 9 the time histories for the microphone signals 
were processed using the finite Fourier transform with suitable averaging and scaling to obtain narrow band 
spectral data with amplitude levels adjusted to report results at a radius of 100 nozzle exit diameters. The 
atmospheric attenuation inherent in the measured data was removed resulting in lossless spectra. These 
spectra contained 8193 points at equally spaced intervals of about 12 Hz. The spectral amplitudes were 
normalized by the frequency interval to obtain a power spectral density. Follow-on processing generated 
1/12-octave and 1/3-octave spectra as desired. 

Given that only the total radiated jet noise was measured, Khavaran & Bridges 9, 10 developed scaling 
laws that attempt to separate the total jet noise into its mixing noise and shock-associated noise components. 
They also discuss the phenomenon of mixing noise amplification resulting from the presence of any shock- 
related screech tones. This additional noise could be considered to be part of the mixing noise component or 
the shock-associated noise component. Dahl 2 showed that the spectra predicted from models for supersonic 
jet mixing noise were closer in comparison to the mixing noise component spectra with amplification noise. 
Hence, the shock-associated noise component spectra used in this study contained no amplification noise. 

The results of the spectral component separation for the Table 1 conditions labeled 2090 and 2095 are 
shown in Figure 3 for measurements at 50, 90, and 120 degree inlet angles. The plots show 1/3-octave 
lossless spectra as a function of the Strouhal number fD/Vj , where, in this case, / is the 1/3-octave center 
frequency, D is the nozzle exit diameter, and V :] is the fully expanded jet velocity. Mixing noise can be seen 
to dominate at the lower frequencies and become a larger portion of the spectrum as the inlet angle moves 
toward downstream. Conversely, the broadband shock-associated noise dominates at the higher frequencies 
and at the upstream angles. The process of separating jet noise spectral components involves using scaling 
laws and interpolations that contain some inaccuracies. One example is shown in Figure 3f where shock 
noise is low relative to the total noise at the downstream angle. The separation method inaccuracies result 
in the mixing noise spectrum having some levels higher than the total noise spectrum levels. Figure 3 also 
shows the predicted spectra from the SAE broadband shock-associated noise model that is given next. 

IV. Description of Prediction Model 

Supersonic jet flows from convergent nozzles are imperfectly expanded, creating shocks in the flow. The 
turbulence normally created by jet flows interacts with the shock structure producing a noise source in 
addition to that produced by turbulent mixing. This source may emit tones, known as screech, and a wide 
band of noise called broadband shock-associated noise. This latter noise is of main interest here. Though 
considered a broadband noise, shock-associated noise has a characteristic spectrum with a fundamental peak 
and has spectral levels that vary with the inlet angle. It is also most dominant at angles upstream of the nozzle 
exit. A prediction model for shock-associated noise was first proposed by Harper-Bourne & Fisher . 12 The 
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1/3rd-0ctave SPL (dB) - Lossless 1/3rd-Octave SPL (dB) - Lossless 1/3rd-Octave SPL (dB) - Lossless 



(a) Mj = 1.392, Tt/Too = 1.76,/? = 0.968,6 = 50 



(c) Mj = 1.392, Tt/Too = 1.76,0 = 0.968,6 = 90 



(b) Mj = 1.659, Tt/Too = 2.16,/? = 1.324,6 = 50 



(d) Mj = 1.659, Tt/Too = 2.16,/? = 1.324,6 = 90 



(e) Mj = 1.392 , Tt/Too = 1.76,0 = 0.968,6 = 120 



(f) Mj = 1.659, Tt/Too. = 2.16,/? = 1.324,6 = 120 


Figure 3. 1/3-Octave spectra of the total, mixing component, and shock component noise for measured data from 
two jets with shocks at 50, 90, and 120 degree inlet angles and distance 100H. The predictions from the original SAE 
broadband shock-associated noise model are shown for comparison. 
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model is discussed further and compared to measured data in Balsa et al. 13 and subsequently documented 
in SAE ARP 876E. 3 The model was eventually incorporated into NASA’s Aircraft Noise Prediction Program 
(ANOPP) 4 as the SAESHK module. After the equation for the far field pressure spectrum is integrated to 
obtain the sound pressure spectrum in 1/3-octave bands, as described by Balsa et ah, 13 the equation for the 
far-held mean-square acoustic pressure ( p 2 ) in each band can be expressed as 


<P 2 > 


= 1.92 x 10" 3 


D 2 ^(PooCgo ) 2 
16 r 2 [1 — Moo cost 1 ] 4 


H(a) [1 + W(a,e,Vj)} 


where the parts of the equation are described as follows: 


( 6 ) 


D - exit diameter of a converging nozzle in meters, 
r - distance from nozzle exit to observer in meters. 

/? - shock noise parameter defined in equation (4). 

ij - exponent that is equal to 4 if (3 < 1 or if fj > 1, it equals 1 for T t /T ^ < 1.1 or 2 for 

T t /Too >1.1 

Poo - the ambient density in kg/m 3 . 

Mqo - the ambient Mach number. 

9 - polar angle starting from the inlet direction. 
a frequency parameter given by 


fD 

a = 2ttK 0 - — /3[ 1 — Moo cos0] 

C-OC) 


( 7 ) 


/ frequency in hertz. 

H - group source strength spectrum. For cold jets T t /T, ' x < 1.1, use log (H) — 0.2. 
W - shock cell interference function given by 


. JV.-l N s -(k+l) . 

jrr 4 V-^ rrx \1 k 2 ST' Sin {bcrqkm/2) I X 

W = WhZ-^ [ A 2 C ( f7 )] zL T7 C0S Mfcm) 


k = 1 


m— 0 


Qkm 


Ah 


KoCoJ (Vj/Coo) 


1 — e l m + 


k + 1 


1 + On — cos 


N s - number of shocks, variable, typically set to 8. 
b - 1/3-octave bandwidth constant, 0.23077. 

C - correlation coefficient spectrum. 

I< 0 - average shock cell spacing, set to 1.1. 

Ah - first shock cell spacing constant, set to 1.31. 

AT 2 - constant multiplying the correlation coefficient spectrum, set to 1.0. 
Cq - eddy convection velocity coefficient, set to 0.7. 
e - due to variation in shock spacing, set to 0.06. 


(8) 

( 9 ) 


V. Results for SAE Model Predictions 

As mentioned above, Figure 3 contains the corresponding spectrum from the SAE model designed to 
predict broadband shock-associated noise. As can be seen in the comparisons at 50 degrees in Figure 3(a), 
the SAE model predicts the shock-associated noise quite well for the lower-speed jet. At the higher speed, 
Figure 3(b) shows the predicted spectrum compares favorably with the total noise spectrum especially at 
lower frequencies where mixing noise is dominant, but not quite as well at frequencies higher than the 
peak frequency where shock-associated noise dominates. At 90 degrees in Figure 3(c) and 3(d), the higher 
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frequency predictions for the shock-associated noise are lower creating a larger difference between predicted 
and measured spectra. In addition, a low frequency hump develops in the predicted spectra that does not 
exist in the shock noise component spectra and does not represent the total noise spectral shape, either. 
This low frequency hump further develops at 120 degrees, as shown in Figure 3(e) and 3(f), and comes 
to dominate the spectra over that of the shock-associated noise related peak as /? and the jet temperature 
increase. There is little agreement between the SAE model predictions and the measured data at these 
conditions. 

The major parts of the SAE model equation ( 6 ) are plotted in Figure 4. The plots show the shape of the 
H- spectrum multiplied by an amplitude coefficient when IT = 0 is used in equation ( 6 ). This spectrum is 
not a function of the inlet angle with no ambient flow and the peak shifts to lower frequencies as /3 increases 
which follows from equation (7). The interference function W modifies the shape of the if -spectrum. Since 
W goes to zero at high frequencies, for plotting purposes the logarithm of 1 + IT added to 100 is displayed 
in the figure. The purpose of the interference function, in effect, is to suppress the group source strength 
ii-spectrum at the lower frequencies and to accentuate the shock-noise peak. Clearly as the inlet angle 
increases and the jet becomes hotter, the group source strength spectrum becomes dominant and W has less 
of an impact. The shock-associated noise is no longer predicted by the model. 

VI. Toward a Revised Model 

A study was performed to seek improvements in the SAE broadband shock-associated noise prediction 
model using the data set given by Table 1 and the ability to separate the measured spectral data into 
component parts related to mixing noise and shock-associated noise. Tester et al . 14 performed a similar 
type of study of this shock-associated noise model. They performed a parametric study on the frequency 
scaling, the correlation coefficient spectrum, the number of shocks, the eddy convection velocity constant, 
and the shock spacing constants. The conclusions based on comparing predicted spectra with total jet noise 
spectra where shock-associated noise dominates were to confirm the choice of parameter values given by 
Harper-Bourne & Fisher 12 and listed following equation ( 6 ) above. Modifications to the model were made in 
computing the overall levels where shock noise levels were extracted from the total noise and differences were 
noted between noise levels emitted by cold jets and levels emitted by jets that were heated. The approach 
here was to study the overall levels using the current approach for separating the shock-associated noise 
component from the measured total noise, conduct a sensitivity analysis of the model parameters, and to 
perform a least-squares fit of the model to the shock-associated noise component spectra of the measured 
data. 

VI. A. Overall Levels 

Equation ( 6 ) shows the results from the model development that the overall levels of shock-associated noise 
vary as (3 to some power 77 . Figure 5 shows plots of the overall sound pressure levels (OASPL) determined 
for all the jet cases shown in Table 1 grouped by jet total temperature ratio. The data are the OASPL at 
three inlet angles less than or equal to 90 degrees. Figures 5(a) to 5(d) use total noise OASPL data. In all 
cases for (3 less than about 1.0, the data follow a f3 4 sloped line except for the lowest (3 cases of the heated 
jets where mixing noise becomes important at these upstream angles. For f3 greater than 1.0, the cold jet 
data in Figure 5(a) follow a /3 1 line and the hot jet data in Figures 5(b) to 5(d) follow a (3 2 line. These result 
are the same as those obtained by Tester et al . 14 where shock noise levels separated from the total noise 
levels were used to perform the scaling. 

The results for the OASPL scaling using the shock-associated noise component determined using the 
Khavaran & Bridges 10 method of component separation is shown in Figures 5(e) to 5(h). The cold jet 
results in Figure 5(e) are the same as in Figure 5(a) for the total noise OASPL. However for the heated 
jet OASPL determined from only the shock-associated noise component, the scaling results for 77 deviate 
from those derived from using the total noise OASPL. At this point, it is possible to best fit a line through 
the data points that would result in non-integer values for 77 . Alternatively, we could follow Khavaran & 
Bridges 10 and allow 77 to be a function of temperature and directivity angle. Most likely, the model is not 
that accurate that these small changes in 77 would be significant in any particular prediction. Hence, we 
use integer values for 77 to maintain the main variation in OASPL with /3, as followed by Harper-Bourne & 
Fisher , 12 Balsa et al ., 13 and Tester et al . 14 For (3 less than about 1.0, the OASPL data follow a /3 5 sloped 
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1/3rd-0ctave SPL (dB) - Lossless 1/3rd-Octave SPL (dB) - Lossless 1/3rd-Octave SPL (dB) - Lossless 



Strouhal Number 


(a) Mj = 1.392, Tt/Too = 1.76, f3 = 0.968, <9 = 50 



(b) Mj = 1.659, Tt/T^ = 2.16, (3 = 1.324 , 6 = 50 



Strouhal Number 


(c) Mj = 1.392, Tt/Too = 1.76, /3 = 0.968, 6 = 90 


Q. 

(/) 


60 4 ; 

10 


10 10 
Strouhal Number 


10 


(d) Mj = 1.65 9, Tt/Too = 2.16,0 = 1.324 ,6 = 90 





(e) Mj = 1.392, Tt/Too = 1.76, 0 = 0.968,0 = 120 (f) Mj = 1.659, Tt/T*, = 2.16,0 = 1.324,0 = 120 


Figure 4. 1/3-Octave spectra of the SAE shock noise prediction and its component parts compared to the shock noise 
component of measured data at 50, 90, and 120 degree inlet angles and distance 1007). 
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(a) Total Noise, Tt/Too = 1.01 (b) Total Noise, Tt/Too = 1.77 




(c) Total Noise, Tt/Too = 2.16 (d) Total Noise, Tt/Too = 2.66 




(e) Shock Noise, Tt/Too = 1.01 (f) Shock Noise, Tt/Too = 1-77 




(g) Shock Noise, Tt/Too = 2.16 


(h) Shock Noise, Tt/Too = 2.66 


Figure 5. Variations of overall sound pressure levels (OASPL) for the jets with operating conditions listed in Table 1 
at three inlet angles as a function of the shock parameter /3. Comparison of OASPL for the jet total noise to OASPL 
for the jet shock noise component only. Solid lines for variation in OASPL versus /3 71 in the SAE model. Dashed lines 
are revised variations in shock noise component only OASPL versus ( 3 11 . 


N ASA/TM— 20 14-218109 


9 


line and for /3 greater than 1.0, hot jet OASPL follow a /3 1 line, the same as the cold jet data. The behavior 
of the shock noise only OASPL at 50 degrees and /3 > 1 was noted but not studied further here. 

VI. B. Sensitivity Analysis 

A sensitivity analysis was conducted of the parameters that affect the behavior of the shock cell interference 
function W. Specifically, the following parameters were varied over the given values: 

A'o = 0.6, 0.9, [1.1], 1.2, 1.3 
A'i = 1.01, 1.21, [1.31], 1.41, 1.61 

e = 0.04, 0.05, [0.06], 0.07, 0.08 (10) 

K 2 = 0.8, 0.9, [1.0], 1.1, 1.2 
C 0 = 0.5, 0.6, [0.7], 0.8, 0.9 

where the terms in brackets are the original model parameter values. Each parameter was varied one at a 
time while the others were held constant at their original value. The effects on 1 + IP are shown in Figure 6 
for the jet condition labeled 2097 where the red dashed lines are the results for the minimum parameter 
values given in equation (10). The changes are systematic to the maximum parameter value results shown 
by the thick black dashed line. The measured data and the component parts of the original model as shown 
in Figure 4 are also shown in Figure 6. This helps to indicate what changes are required of the shock cell 
interference function W that when it is included with the scaled H-spectrum (W = 0), the model predicted 
spectrum would be closer to the measured spectrum. The changes in the predicted spectrum as W changes 
are not shown to reduce clutter in the plots. 

The variations in 1 + W due to changes in the shock spacing parameters Kq, K\, and e are shown in 
Figures 6(a) to 6(c). The largest changes are due to K\ with lesser effects due to changes in e and K 0 - 
Note that changing K\ affects the location of the shock-associated noise peak. As discussed in Tester et 
ah, 14 these changes could improve the spectrum prediction at some frequencies, but degrade the prediction 
at other frequencies. Adjusting these parameters to fit predicted spectra to measured data could lead to the 
shock spacing parameters having nonphysical values. 

The K 2 parameter is intended to increase or decrease the correlation coefficient spectrum C(<r) value 
without changing the shape of this spectrum. C(cr) has the highest correlation values at the lower frequencies 
with decreasing values as the frequency increases. Using the documented values of C(ct) as a baseline, 3,4 
K 2 is limited to a maximum value of 1.351 to insure that the correlation coefficient spectrum is less than 
or equal to 1.0. The results in Figure 6(d) shows that increasing K 2 would result in a slightly more rapid 
decrease in the lower frequency portion of the predicted spectrum. 

Finally, Figure 6(e) shows that changes in the eddy convection velocity coefficient Co has the largest 
effect on 1 + W . It will be considered further separately. 

VI. C. Convection Velocity 

The results shown in Figure 6(e) indicate that changes in the convection velocity affect the behavior of 1 + W 
in the prediction model. These changes are for a constant convection velocity at all frequencies. As with the 
shock spacing parameters, its adjustment could improve or degrade the predicted spectrum depending on the 
frequency. But unlike the shock spacing parameters which are fixed values that should be determined from 
flow measurements, the eddy convection velocity coefficient is known to be frequency and velocity dependent 
through the Strouhal number. 15-17 

To explore the effect of a Strouhal dependent convection velocity on the predicted shock-associated noise 
spectrum, the Co eddy convection velocity coefficient in equation (9) was considered to be a function of 
frequency. A least-squares fit of the predicted spectrum to the measured shock-associated noise component 
spectrum was performed where the Co’s at each frequency are the only parameters for each of 15 hot jet cases 
where My > 1.05 and where the fit included data from eight inlet angles from 50 to 120 degrees. Using the 
velocity V) for each operating condition, the results for Co are plotted as a function of the Strouhal number 
in Figure 7. The Strouhal number range was limited to about 0.2 to 2.0 to cover the peak shock-associated 
noise. Even though the fit is done using far-field spectral data, the eddy convection velocity coefficient Co 
values tend to cluster along the lines given by Harper-Bourne 15 and Morris & Zaman 16 for the Strouhal 
dependent convection velocity derived from flow measurements especially for Strouhal numbers less than 1.0 
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1/3rd-0ctave SPL (dB) - Lossless 1/3rd-Octave SPL (dB) - Lossless 




(a) Vary Kq, average shock cell spacing (b) Vary K i, first shock cell spacing constant 



(c) Vary e, variation in shock spacing 



(d) Vary K2, correlation coefficient spectrum multiplier 



(e) Vary Co, eddy convection velocity coefficient 


Figure 6. 1/3-Octave spectra results from a sensitivity analysis of the parameters affecting the 1 + W component of the 
SAE shock noise model. Jet condition 2097 at a 90 degree inlet angle and distance 100D. 
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Figure 7. Scatter of convected velocity results from a least-squares fit of the SAE shock noise model using the shock 
noise component of measured far-field noise data. Lines are the results from measured flow field data. 


and P > 0.967. In both of the previous measurement cases, low speed jets were used. Dahl 17 later showed 
convection velocity results from a near M ? = 1 jet that followed these convection velocity lines. Hence, we 
replace the constant Co by the equation from Morris & Zaman 


Co (/, Vj) = 0.062 In 



0.701 


( 11 ) 


VI. D. Group Source Strength Spectrum 

Of the parameters available to improve equation (6) to better predict the shock-associated noise component 
of the measured data sets, we are left with the group source strength spectrum H(a) and the correlation 
coefficient spectrum C(a) at each frequency. Attempts at adjusting both in a least-squares fit has not led 
to satisfactory results as the best fit values at each frequency had large variations from one frequency to the 
next and from one operating condition to the next. Hence, the current choice is to set K% = 1.0 and to use 
the original model values for C(a) as suggested by the findings of Tester et al. 14 

Using equation (11) for the eddy convection velocity coefficient and the revised set of r) exponent values 
for hot jets given in Figures 5(f) to 5(h), a least-squares fit of the predicted spectrum to the measurement 
derived shock-associated noise component spectrum was performed to obtain H(a). The data was from the 
20 jet cases shown in Table 1 where M ? - > 1.05 including five cold jet operating conditions and the remainder 
were hot jets. The spectral data came from three inlet angles, 50, 60, and 70 degrees, where shock-associated 
noise dominates at all operating conditions. The fit results for H(a) are shown in Figure 8(a) for all the 
jet operating conditions. The original ff(cr) spectrum is shown for comparison. 3,4 The remaining five plots 
in Figure 8 show the same data at nominally fixed /3 values representative of those shown in Table 1. The 
findings are as follows: 

• The fit results for H(cr) at frequencies less than the peak frequency have larger variations due to 
the shock-associated noise component spectra containing oscillations and possible remnants of screech 
tones as can be seen in the shock noise component spectra in Figure 3. 

• H(a) appears to be independent of jet temperature for hot jets at a fixed /3 as noted by the similar 
results at frequencies higher than the peak frequency. 

• The cold jet results for H(a) were found to be about 3dB lower than the hot jet results for H(a). The 
cold jet H(a) values in Figure 8 are plotted as log(i7) + 0.3. This differs from the current SAE model 
where cold jet levels differ from hot jet levels by 2dB. 3,4 
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log(H) log(H) log(H) 








Figure 8. Results for the group source strength spectrum H from a least-squares fit of the SAE shock noise model 
using the shock noise component of measured far-field noise data. Lines are the original SAE model if-spectrum. 
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• The issue with the oscillation in the measured shock-associated noise component spectra at frequencies 
less than the peak frequency (mentioned above) are especially evident in the /3 = 0.628 jet results, 
Figure 8(a) (red symbols). These lower frequency points are not considered to be representative of the 
lower frequency decrease in shock-associated noise. Removing these points shows that the remaining 
H(a) fit results have a definite clustering about a new spectrum shape. Furthermore, the removing of 
the /3 = 0.628 fit points does not have a significant effect on the peak and higher frequency definition 
of a new FF(cr) spectrum. 

The objective is now to replace the SAE FF-spectrum with a new spectrum based on the fit results for 
FF(er) where /3 > 0.967. This was done using the empirical mode decomposition method and using the 
least-squares fit method. 

VI.D.l. H(a) from Empirical Mode Decomposition 

The cluster of FF(cr) fit data points shown in Figure 8(a) may be considered as a short duration signal with 
high frequency content riding on a definitive background trend. A method to extract this trend is called 
empirical mode decomposition. This method separates a data signal into a series of basis functions each 
containing a limited frequency range starting from a high frequency band and continuing down toward lower 
frequency bands until a residual of the signal is reached that contains no oscillating content. 18 For our 
purposes, the residual of the series of data points put together from the FF(ct) fit points represents the basic 
trend or shape of the FF(cr) spectrum. Dahl 1 ' described the process of empirical mode decomposition (EMD) 
that was used on this data. Given that the FF(cr) signal is made from overlapping spectral data points, 
the signal is highly irregular in its sampling. Since EMD was applied to this set of data points regardless 
of variable sampling rate (Hence, implicitly evenly spaced.), the resulting smooth residual curve contained 
jaggedness when plotted versus the actual log(cr) value for each data point. The solution to obtaining a 
smooth log (FF) versus log(cr) curve was to pass the residual signal through a 15-point moving average filter 
and then resample it onto an evenly spaced grid in log(tr) using a 5-point Lagrange interpolating polynomial. 
To extend the curve beyond the high end of the data out to log(cr) = 2.5, a straight line was used based 
on a linear least-squares fit of the EMD spectrum between log(cr) = 1.6 to 1.72. The process of computing 
the EMD involves choosing stopping criteria that may lead to variations in results. A sampling process was 
employed to reduce the variation by averaging the results. After 200 samples, the averaged FF-spectrum had 
a maximum uncertainty of ±0.005. This averaged EMD-derived FF-spectrum is shown in Figure 9 compared 
to the data points and the original SAE FF-spectrum, and it is also tabulated in Table 5, Appendix B. 



Figure 9. Single curve definitions for the if-spectrum derived from fit data. 
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VI. D. 2. Equation for H(a) from the least-squares fit method 

The H(a) fit data points shown in Figure 9 also have a shape reminiscent of the shock-associated noise 
spectral shape proposed by Deneuville 19 consisting of two sloping straight lines meeting at the peak of the 
spectrum. But instead of an abrupt change in slopes between two lines, it was desired to represent the smooth 
curvature at the peak that EMD showed was part of the residual trend of the spectrum. The solution was 
to connect the two straight lines with two lines governed by cubic equations. At the three interface points 
between the four lines, the value and its derivative were considered to be continuous. It turns out that a 
sensible least-squares fit is difficult to achieve if the location of the three interface points are included as fit 
parameters. Thus, the choice of the location of these three points was arbitrary. 

Using x = log(a) and y = log (H), the 4 line piecewise curve fit to the H(a) fit data is given by: 

x < 0.5 

y = —1.431 + 4.284 {x - 0.5) 

0.5 < x < 0.75 

y = -0.6307 + 0.01952 (x - 0.75) - 21.13 (x - 0.75) 2 - 33.60 (x - 0.75) 3 

(12) 

0.75 < x < 0.9 

y = -0.6307 + 0.01952 (x - 0.75) - 4.052 (x - 0.75) 2 + 5.882 (a: - 0.75) 3 
0.9 < x 

y = -0.6991 - 0.7989 (x - 0.9) 

where the maximum uncertainty in the fit of y is ±0.1. A plot of this spectrum is also shown in Figure 9. This 
spectrum result is very close to that obtained using EMD except at the highest frequencies. The extrapolation 
of the linear fit to frequencies higher than log(cr) = 1.72 will result in higher curve-fit spectrum levels than 
for the EMD-derived spectrum extrapolated levels. 

VII. Revised Model Results 

The predicted spectra from the revised model are compared to the original SAE shock-associated noise 
prediction model. Using the same operating conditions and inlet angles as in Figure 3, the revised model 
shock-associated noise predictions are shown in Figure 10 for both the EMD-derived and the curve fit H - 
spectra. As expected from the closeness in the new 77-spectrum shapes shown in Figure 9, the results in 
Figure 10 are nearly identical except near the extremes of high and low frequencies. The revised predictions 
more closely represent the shock-associated noise component spectra of the measured data at inlet angles 
less than or equal to 90 degrees. The lower frequency hump arising in the original SAE model spectra is no 
longer present. At the 120 degree angle, the revised model is also failing in the same manner as the original 
model even though the revised model is an improvement by lowering the levels at the lower frequencies and 
raising the levels at the higher frequencies. The group source strength spectrum still becomes dominant and 
the interference function is unable to cancel out the lower frequency hump and accentuate the shock noise 
peak sufficiently for better agreement with the measured data. 

One consequence of using the new i7-spectra shown in Figure 9 that are based on measured data from 
jets with higher ft values is that jets with lower ft values may have their shock-associated noise component 
spectra inaccurately predicted. In Figure 8(a), it is shown that the ft = 0.628 spectral data points are higher 
in level at the lower frequencies than the higher ft jet spectral data points. Thus, the revised model will 
underpredict the data at the lower frequencies. (The original SAE model has similar spectral characteristics 
as a function of ft as shown in the large set of comparisons between predicted and measured spectra in 
Dahl. 2 ) However, the lower ft jets have similar levels to those of the higher ft jets at the peak frequency 
and higher. Hence, the higher frequencies should still be well predicted. This fact will be accounted for in 
computing the prediction errors. 

The errors and uncertainties were computed for the original SAE model and the revised model shock- 
associated noise predictions at the 24 jet operating conditions listed in Table 1. A measure of the error in 
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1/3rd-0ctave SPL (dB) - Lossless 1/3rd-Octave SPL (dB) - Lossless 1/3rd-Octave SPL (dB) - Lossless 




(a) Mj = 1.392, Tt/Too = 1.76, 0 = 0.968 , 6 = 50 (b) Mj = 1.659, T t /Too = 2.16, 0 = 1.324, 0 = 50 




(c) M 3 = 1.392, Tt/Too = 1.76, 0 = 0.968, 6 = 90 (d) Mj = 1.659, Tt/Too = 2.16,0 = 1.324, 0 = 90 




(e) Mj = 1.392, Tt/Too = 1.76,0 = 0.968,6/ = 120 (f) Mj = 1.659, Tt/Too = 2.16,0 = 1.324,0 = 120 


Figure 10. 1/3-Octave spectra of the shock-associated noise component for measured data from two jets with shocks 
at 50, 90, and 120 degree inlet angles and distance 100T. The predictions from the original SAE broadband shock- 
associated noise model, the revised model with the EMD-derived 17-spectrum, and the revised model with the curve 
fit 17-spectrum are shown for comparison. 
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the predicted spectrum was computed by averaging the errors computed for each 1 /3-octave band 

1 N 

Average Error = — ^ |SPL pred (/, ; ) - SPL meas (/i)| (13) 

V i—1 

where N is the number of included 1/3-octave bands. By taking the absolute value of the difference between 
predicted and measured sound pressure levels, a smaller average error requires relatively smaller errors across 
the entire spectrum. Another measure of error was computed using the overall sound pressure level (OASPL) 
using 

OASPL Error = OASPL pred - OASPL meas (14) 

The calculation of OASPL for this error determination was done over the same frequency range as that 
used to compute the average error. The OASPL error can be either positive or negative indicating that 
the predicted value is higher or lower, respectively, than the measured OASPL. These errors were computed 
separately for each inlet angle. 

The errors were computed using equations (13) and (14). For the nominally p = 0.314 cases, the error 
computation was limited to the frequency range where the Strouhal number was greater than or equal to 
1.0. This was done because the model underpredicts the data at lower frequencies for lower f3 conditions 
but is still fairly accurate at the peak of the shock-associated noise component spectrum and at higher 
frequencies. The lower frequencies of the total jet noise spectrum are due to mixing noise sources. Hence, 
the sum of the mixing noise model prediction and the shock noise model prediction in lower P cases will give 
a reasonable prediction for the total noise spectrum. The /3 = 0.628 cases were limited to Strouhal numbers 
greater than or equal to 0.4. All other higher p cases had no Strouhal number range limitations. Plots of 
the errors as a function of the inlet angle are shown in Figure 11 for the original SAE model and for the 
revised model using the EMD-derived H-spectrum. The error results for the curve fit essentially directly 
overlap the EMD-derived results and therefore are not shown. The two lower p results are shown separately 
due to their frequency limitations. These results are given in terms of the range of errors computed over the 
four cases at that p value. The remaining 16 cases are plotted together with the values at each inlet angle 
averaged to obtain an average error within the range of computed errors. 

The results for average error and OASPL error show that the revised model predictions generally has 
lower errors than the original SAE model predictions. Both models have typically larger errors at inlet angles 
greater than 100 degrees. The data in Figure 11 was used to determine an average decrease in error from 
the SAE model to the revised model. Also, the range of errors was used to compute the uncertainty in both 
model predictions. These results are shown in Table 2. The table shows that the revised model prediction 
is as good as the SAE model prediction with zero decrease in error to up to a 3.5 dB decrease in the average 
error for the revised model. For the most part, the revised model predictions have less uncertainty than the 
SAE model predictions. 

Table 2. Average of errors and uncertainties computed from comparisons of predicted to measured shock-associated 
noise component spectra for both the SAE model and the revised model. 



Average Error 

OASPL Error 



Avg. Err. 

SAE Mod. 

Rev. Mod. 

Avg. Err. 

SAE Mod. 

Rev. Mod. 



Decrease 

Uncert. 

Uncert. 

Decrease 

Uncert. 

Uncert. 

Conditions 

P = 0.314 

0.0 

±1.04 

±0.64 

2.0 

±4.2 

±3.3 

e < ioo, fD/Vj > i.o 

P = 0.628 

0.5 

±0.62 

±0.33 

0.7 

±1.4 

±1.4 

e < 120, fD/Vj > 0.4 

P > 0.967 

3.5 

±1.73 

±1.36 

1.4 

±2.1 

±1.7 

e < 120, fD/Vj > 0.0 


Errors and uncertainties in decibels. 


The error results listed above are for the data used to develop the revised shock-associated noise com- 
ponent prediction model. The revised model is now applied to jet conditions not included in the database 
of Table 1. The two supersonic jet conditions are listed in Table 3; a heated jet at T t /T, x = 1.45, less than 
the heated jet values in the database, and a higher temperature jet with Tt/T, = 3.20 using the smooth 
srncOOO nozzle with low screech. The noise predictions were made with the inclusion of mixing noise from 
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Average Error in dB Average Error in dB Average Error in dB 



Inlet Angle 



(a) p = 0.314 


(b) (3 = 0.314 



Inlet Angle 



Inlet Angle 


(c) p = 0.628 


(d) p = 0.628 



Inlet Angle 


(e) p > 0.967 



Figure 11. Summary of error between model predictions and measured data. Average error across the spectrum for 
the shock noise component on the left. Error in OASPL on the right. Results from prediction models minus measured 
data. 
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Table 3. Test conditions for a convergent nozzle not listed in Table 1. 


Nozzle 

Tt/T,,, 

NPR 

Mi 

Vj / ^OO 

P 

7 j 

Too(K) 

Poo (Pa) 

Coo (m/s) 

smc021 

1.45 

3.653 

1.497 

1.500 

1.114 

1.3997 

274.9 

99811 

332.31 

srncOOO 

3.20 

4.017 

1.568 

2.300 

1.208 

1.3752 

276.3 

99794 

333.17 


the prediction model in SAE ARP 876E. 3 This shows how well the models predict the total noise spec- 
trum where shock-associated noise dominates at the upstream angles and whether or not the shock noise 
prediction model’s lack of accuracy matters as the angle moves toward downstream where mixing noise is 
more important. The comparison of the total noise prediction to measured data is shown in Figure 12 for 
both the SAE shock noise model prediction with the added mixing noise prediction from the SAE mixing 
noise model and the revised shock noise model prediction also with the added mixing noise prediction. At 
90 degrees and upstream, the revised model clearly shows better spectral predictions, limiting the effect of 
the low frequency hump evident in the SAE model predicted spectrum and increasing the spectral levels at 
the peak and higher frequencies. For the cooler jet at the 120 degree inlet angle, Figure 12(e) shows the 
revised model gives a slightly better prediction while Figure 12(f) shows the hot jet spectrum is dominated 
by mixing noise and the accuracy of the shock-associated noise prediction is not relevant. The computed 
errors are shown in Table 4. The errors are lower for the revised model compared to the SAE model where 
shock-associated noise dominates at the upstream angles and about the same where mixing noise dominates. 

Table 4. Computed errors from comparison of predicted to measured total jet noise spectra shown in Figure 12. 



Average Error 


P = 

1.114 

P = 

1.208 

Inlet 9 (Deg) 

SAE Model 

Rev. Model 

SAE Model 

Rev. Model 

50 

4.3 

2.8 

2.8 

1.5 

90 

4.3 

2.5 

3.3 

2.0 

120 

4.0 

3.0 

2.8 

2.6 


OASPL Error 

50 

-1.90 

-0.10 

-1.52 

-0.07 

90 

-2.73 

-0.62 

-2.37 

-1.55 

120 

-0.01 

+0.43 

+2.19 

+2.15 


Errors in decibels. 


In summary, the proposed changes to the SAE shock-associated noise prediction model expressed in 
equations (6) to (9) with the associated parameter definitions are as follows: 

1. The /? exponent r\ is equal to 1 if /3 > 1 or if fi < 1, it equals 4 for Tt/Too < 1.1 or 5 for T t /T, > 1.1. 

2. The eddy convection velocity coefficient Co is given as a function of the Strouhal number by equation 

( 11 ). 

3. The group source strength spectrum H is given in terms of log(H) versus log(cr) by either interpolation 
from the values in Table 5, Appendix B, or by the curve fit equations (12). The tabulated values come 
from the residual of an empirical mode decomposition that gives the basic trend or shape of the 
spectrum with no a priori assumption about the spectrum shape. The piecewise curve fit equations 
and their points of connection were arbitrarily chosen to fit a curve to represent the data. 

4. The cold jet results are 3dB lower than the hot jet results. Replace tabulated or computed log(H) 
with log (H) — 0.3 when T t /T : < 1.1. 
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(a) Mj = 1.497, Tt/Too = 1.45, jS = 1.114,0 = 50 



(c) Mj = 1.497, Tt/Too = 1.45, /} = 1.114,0 = 90 



(e) Mj = 1.497, Tt/T^ = 1.45, ,3 = 1.114,0 = 120 



(b) Mj = 1.568, Tt/Too = 3.20,/J = 1.208,0 = 50 



(d) Mj = 1.568, Tt/Too = 3.20 , 0 = 1.208,0 = 90 



Strouhal Number 


(f) Mj = 1.568, Tt/Too = 3.20, 0 = 1.208,0 = 120 


Figure 12. 1/3-Octave spectra of the total noise for measured data from two jets with shocks at 50, 90, and 120 degree 
inlet angles and distance 1000. The predictions from the original SAE broadband shock-associated noise model and 
the revised model with the EMD-derived F-spectrum are shown for comparison. Both models are summed with the 
SAE mixing noise model to obtain the total predicted jet noise spectra. 
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Appendix A — Computation of Ratio of Specific Heats 


The ratio of specific heats in a supersonic jet at the fully expanded jet conditions depends on the jet 
temperature. Assuming that the jet flow is a calorically imperfect, thermally perfect gas, it is computed by 
iterating on the following equations: 20 


My(j) — 


L7j(i— 1) 


- 1 


NPR "W-d - 1 


(15) 




T t 


7j(i- 1) 


-n-1 




7j(») - 1 + 


Tperf 1 

1 (Tperf 1) 

/ 0 \ 2 e e /7j(0 

\ T m) - 1) 2 _ 


(16) 


(17) 


where the nozzle pressure ratio, NPR, and total temperature, T t , are known and fixed operating conditions, 
7j(0) = 7 P erf = 1-4 for air treated as a perfect gas, and 


0 = 


5500. , T in degrees Rankine 

5500. /1. 8, T in degrees Kelvin 


(18) 
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Appendix B — //-Spectrum table 


Table 5. Revised group source strength spectrum H(c r). 


logio a 

logio H 

logio G 

logio H 

b 

o 

bJD 

JO 

logi 0 H 

b 

o 

bio 

JD 

logio H 

0.000 

-3.5655 

0.625 

-0.8350 

1.250 

-0.9419 

1.875 

-1.5835 

0.025 

-3.4736 

0.650 

-0.7547 

1.275 

-0.9578 

1.900 

-1.6118 

0.050 

-3.3813 

0.675 

-0.7028 

1.300 

-0.9763 

1.925 

-1.6400 

0.075 

-3.2648 

0.700 

-0.6697 

1.325 

-0.9989 

1.950 

-1.6683 

0.100 

-3.1415 

0.725 

-0.6528 

1.350 

-1.0231 

1.975 

-1.6966 

0.125 

-3.0261 

0.750 

-0.6488 

1.375 

-1.0477 

2.000 

-1.7249 

0.150 

-2.9135 

0.775 

-0.6531 

1.400 

-1.0742 

2.025 

-1.7532 

0.175 

-2.8252 

0.800 

-0.6638 

1.425 

-1.0978 

2.050 

-1.7815 

0.200 

-2.7293 

0.825 

-0.6776 

1.450 

-1.1205 

2.075 

-1.8098 

0.225 

-2.6481 

0.850 

-0.6952 

1.475 

-1.1399 

2.100 

-1.8381 

0.250 

-2.5495 

0.875 

-0.7121 

1.500 

-1.1617 

2.125 

-1.8663 

0.275 

-2.4462 

0.900 

-0.7323 

1.525 

-1.1841 

2.150 

-1.8946 

0.300 

-2.3262 

0.925 

-0.7519 

1.550 

-1.2107 

2.175 

-1.9229 

0.325 

-2.1857 

0.950 

-0.7731 

1.575 

-1.2366 

2.200 

-1.9512 

0.350 

-2.0466 

0.975 

-0.7921 

1.600 

-1.2677 

2.225 

-1.9795 

0.375 

-1.9187 

1.000 

-0.8103 

1.625 

-1.2991 

2.250 

-2.0078 

0.400 

-1.7948 

1.025 

-0.8271 

1.650 

-1.3346 

2.275 

-2.0361 

0.425 

-1.6939 

1.050 

-0.8417 

1.675 

-1.3605 

2.300 

-2.0644 

0.450 

-1.6027 

1.075 

-0.8542 

1.700 

-1.3857 

2.325 

-2.0926 

0.475 

-1.5233 

1.100 

-0.8667 

1.725 

-1.4074 

2.350 

-2.1209 

0.500 

-1.4252 

1.125 

-0.8786 

1.750 

-1.4420 

2.375 

-2.1492 

0.525 

-1.3137 

1.150 

-0.8910 

1.775 

-1.4746 

2.400 

-2.1775 

0.550 

-1.1790 

1.175 

-0.9017 

1.800 

-1.4986 

2.425 

-2.2058 

0.575 

-1.0566 

1.200 

-0.9146 

1.825 

-1.5265 

2.450 

-2.2341 

0.600 

-0.9327 

1.225 

-0.9267 

1.850 

-1.5552 

2.475 

2.500 

-2.2624 

-2.2907 
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